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SUMMARY 

The equations of motion for large structures with rotating flexible com- 
ponents are derived by regarding the structure as an assemblage of substruc- 
tures. Based on a stationarity principle for rotating structures, it is shown 
that each continuous or discrete substructure can be simulated by a suitable 
set of admissible functions or admissible vectors. This substructure synthe- 
sis approach provides a rational basis for truncating the number of degrees of 
freedom both of each substructure and of the assembled structure. 

INTRODUCTION 

The methodology for analyzing large complex structures has developed along 
different lines. One approach represents a natural extension of methods de- 
veloped originally for civil and aircraft structures, culminating in the fi- 
nite-element method (ref. 1) and the component-mode synthesis (refs. 2,3). Al- 
though rotation of the structure could be accounted for through rigid-body 
modes, work using the approach of references 1-3 has been concerned mainly 
with nonspinning structures. On the other hand, an entirely different approach 
was developed in conjunction with spinning and nonspinning spacecraft struc- 
tures. This approach was dominated by the fact that early spacecraft could be 
treated as entirely rigid. Hence, in the early stages of development, struc- 
tures were assumed to consist of point-connected rigid bodies arranged in "to- 
pological trees" (refs. 4,5). With time, the rigidity assumption was relaxed 
gradually by first allowing for flexible "terminal bodies" (refs. 6,7) and 
then finally for all flexible bodies (ref. 8). A third approach to the prob- 
lem of spinning flexible spacecraft was concerned with spacecraft consisting 
of a rigid body with flexible appendages (ref. 9,10). This latter approach 
can be regarded as an early application of the component-mode synthesis to 
spinning structures. 

Most papers concerned with structures simulated by point-connected rigid 
bodies, such as references 4,5, proposed to derive the equations of motion by 
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the Newtonian approach, on the assumption that such a derivation was more suit- 
able for digital computation. Of course, an early difficulty became immediate- 
ly apparent in the form of the handling of interbody constraints, a major crit- 
icism of the Newtonian approach in most circumstances. Another difficulty was 
the relatively large number of degrees of freedom involved, a difficulty only 
compounded by permitting various bodies to be flexible. As a result, there 
are no meaningful ways of truncating the problem. 

This paper is concerned with the mathematical simulation of large struc- 
tures, where the structure is regarded as an assemblage of substructures. In- 
deed, the mathematical model is assumed to consist of a central substructure 
with a number of appended substructures, where some of the latter can rotate 
relative to the central substructure. To ensure that the various substructures 
act as parts of a whole structure, an orderly kinematical procedure is used 
which takes into account automatically the superposition of motion of the cen- 
tral substructure on the motion of the interconnected substructures. The sys- 
tem equations of motion are derived by means of the Lagrangian approach, which, 
when used in conjunction with the kinematical procedure just described, does 
away with the question of constraints. The equations of motion are derived 
from scalar functions, namely, the kinetic and potential energy, where the 
first requires the calculation of velocities only. In addition, discretiza- 
tion of the kinetic and potential energy in conjunction with linearization 
ensures proper symmetry and skew symmetry of the coefficient matrices in the 
final equations of motion. Using a Rayleigh-Ritz approach, the motion of each 
continuous (discrete) substructure can be represented by a linear combination 
of admissible functions (vectors) rather than substructure natural modes. 

This approach is based on a stationarity principle for rotating structures de- 
veloped recently by the first author (ref. 11). Finally, the truncation prob- 
lem can be handled much more efficiently by the substructure synthesis ap- 
proach, as the possibility of truncating the number of degrees of freedom both 
of the individual substructures and of the assembled structure provides a much 
more rational basis for an overall truncation decision. 

KINEMATICAL CONSIDERATIONS 


Let us consider a general structure consisting of a central substructure 
C and a given number of appended substructures (see fig. 1), where the latter 
are of three types: rigid and rotating relative to the central substructure 

(type R), elastic and nonrotating relative to the central substructure (type 
E), and elastic and rotating relative to the central substructure (type A). 
Clearly, there can be more than one appendage of a given type, but we shall 
confine our discussion to a representative one of each type, with summation im- 
plied over the entire number of substructures. Although we consider here only 
peripheral substructures, the formulation can be easily extended to chains of 
substructures, as discussed later. 

Let us introduce the inertial system of axes XYZ with the origin at 0 and 
identify a system of axes xryczc origin at an arbitrary point C of 

the central substructure. Then, denoting by wqp the radius vector from 0 to C, 
by rQ the position vector of any mass point in the substructure, and by Uq the 
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elastic displacement of that point measured relative to XQyQZ^, and recognizing 
that Wqq is in terms of components along XYZ and r^ and u^ are in terms of com- 
ponents along x^yrZQ, the absolute position of the mass point in question in 
terms of components along x^y^Zj- is Wq = TqqWqc + Tc yc’ where Tqq is the ma- 
trix of direction cosines between XYZ and XQyQZQ. 'Moreover, if fir is the angu- 
lar velocity of the frame XQyQZQ relative to XYZ, the absolute velocity of the 
mass point is 

' ^OC^OC “ 


where Vq + uq is a skew symmetric matrix associated with r^ + Uq and uq is the 
elastic velocity of the point relative to axes XQyQZQ. 


To calculate the absolute velocity of a point in the substructure R, we 
must first obtain the velocity of point R as well as the angular velocity of 
a reference frame XQRy^RZQR attached to the central substructure at R and with 
axes parallel to the rotor axes x^y^zp when at rest and when the central sub- 
structure is undeformed. Due to geometry alone the orientation of axes XQ^y^p 
Zrp relative to XryrZp is given by the constant matrix of direction cosines 
Lqr- Denoting by^ycR^the elastic deformation vector at point R of the central 
substructure and assuming that the components UQp)<_, ucRy» UQpz,, of ucr are small, 
the rotation vector of axes XQRycRZCR due to elastic deformation can be writ- 
ten in the form 


^r^'-grHcr^ 


^^CRz _ ^^CRy 


^^CRx 

^^CR 


du 


CRz 


3x 


CR 


^^CRy _ ^^CRx 
^^CR ^^CR 


( 2 ) 


where is a skew symmetric differential operator matrix corresponding to the 
curl operator. Hence, the matrix of direction cosines between axes x^py^pz^p 
before and after deformation is 

'^^CRx ^\Rz\ 

.'^CR " '^CR / 

( 3 ) 

1 

Moreover, letting Lp be the matrix of direction cosines between axes xpypzp 
and XpnypnZrp, the transformation matrix between axes xpypzp and xcyc^C 
simply Tcr = LpLcpLQp. 

Denoting by top the angular velocity of the rotor relative to axes x^pyQp 
ZQp, the absolute'angular velocity of xpypzp in terms of components along 
xrXrzr is 

~R ^ ^CR ~C ’■R \r^'-GR ^CR^ ~R 
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where the second term in equation (4) is the angular’ velocity of axes XQfjyrRZQR 
due to the elastic motion of the central substructure. Because the rotor is 
rigid, the position of a mass point relative to R is simply r^. Hence, the 
absolute velocity of the point in question is simply 

~R "" ^CR ~CR " *^R ~R 

where Wq|^ is the velocity of point R obtained from Wq by substituting the coor- 
dinates of the point R for those of an arbitrary point. 


Next, let us turn our attention to the substructure E and denote by x^y^zj 
any convenient set of axes with the origin at E and attached to the substruc- 
ture. Using the analogy with equation (4), the angular velocity of x^y^ZE is 


~E ^CE ~C ^CE^^GE ^CE^ 


( 6 ) 


where Tqc = LceLqe* Moreover, by analogy with equation (5), the absolute ve- 
locity of a mass point in the substructure is 

' ■^CE %E - <^E “E>?E ^E 

where u^; is the elastic displacement relative to axes XEyEZE* 

The extension to elastic substructures rotating relative to the central 
substructure is quite obvious. Letting be the angular velocity of the sub- 
structure A relative to a set of axes XQ^yQAZQy\ attached to the central body 
at point A, the absolute angular velocity of x^y^z^ is simply 

~A " ^CA ~C ^ '-A ^CA^4a ^CA^ ^A 

where = L/\Lq/\Lq/\, and the absolute velocity of an arbitrary point in A is 
'^^A " "^CA '5?CA “ ^*"A ^A^?A ^A 


Finally, let us consider chains of substructures. First, we note that 
the angular velocity of a peripheral substructure and the absolute velocity of 
an arbitrary point in a peripheral substructure are written in terms of the 
angular velocity of a set of axes attached to the central substructure and 
with origin at the interconnecting point and the translational velocity of the 
interconnecting point. As an example, see equations (4) and (5). To write 
the angular velocity and absolute velocity of an arbitrary point of a sub- 
structure in a chain, we simply replace quantities pertaining to the central 
substructures, such as Tq[^, ng, Vqr{Lqr ucr)» wqr equations (4) and (5) 
by analogous quantities pertaining to the immediately preceding substructure 
in the chain. 
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SYSTEM DISCRETIZATION AND/OR TRUNCATION 


In general, each elastic substructure possesses a large number of degrees 
of freedom. In fact, if the substructure is continuous, then its number of 
degrees of freedom is infinite. For practical reasons, we must limit the for- 
mulation not only to a finite number of degrees of freedom but also to as small 
a number as possible consistent with a good simulation of the system dynamic 
characteristics. In this regard, we wish to use a Rayleigh-Ritz approach and 
represent the elastic displacements of a continuous substructure by a linear 
combination of space-dependent admissible functions multiplied by time-depen- 
dent generalized coordinates of the substructure. If the substructure is dis- 
crete, then instead of admissible functions we must use admissible vectors. 

Note that it is common practice to use as admissible functions and admissible 
vectors the eigenfunctions and eigenvectors of the substructure. In view of 
the stationarity principle for gyroscopic systems developed in reference 11, 
however, this is not really necessary, and a reasonable set of admissible func- 
tions or admissible vectors should suffice. Hence, we shall use the discreti- 
zation and/or truncation scheme 

yc \ Dc ’ ^E ®E 2e ’ ^A ^ ^A ^A 

where tjq , n^, and n/\ are time-dependent vectors of generalized displacements 
with dimeniions and respectively, and fp, and are 3 x Dq. 

3 X n^, and 3 x Ha space-dependent matrices admissible functions or admis- 
sible vectors, as^the case may be. Note that for a continuous substructure u 
depends on continuous space variables and for a discrete substructure it de-~ 
pends on discrete space variables. In the latter case, the partial derivatives 
involved in the quantity vu are to be replaced by corresponding slopes. 

Although we have mentioned both continuous and discrete substructures in 
the above, we have made no attempt to make clear distinction between the two 
types of mathematical models. Neither have we elaborated on the various types 
3f discrete models, such as lumped models, finite-element models, etc. Of 
:ourse, the mathematical model used depends on the substructure mass and stiff- 
less distributions, but this is of no particular concern here. The reason for 
this is that, independently of the mathematical model postulated for the sub- 
structure, the general idea is the same, namely, to eliminate the spatial de- 
jendence by the use of admissible functions or admissible vectors and to trun- 
:ate the problem by limiting the number of these functions or vectors. 

LAGRANGE'S EQUATIONS OF MOTION 


To derive Lagrange's equations of motion it is necessary to produce first 
ixpressions for the kinetic energy, potential energy, and nonconservative vir- * 
ual work. Assuming that in equilibrium the central substructure C, substruc- 
ure R, and substructure A rotate with the uniform angular velocities fir about 
Q, fi|^ about Zf^, and fi^ about^z^, respectively, while any other motion is zero, 
e can write fi^ " ^ ®C ®C’ ~R ~ ^R ~R ®R §R» ¥A ” ^A ^A ®A ®A’ where 

^ is the vector of direction~cosines between zq and XyZ, Ar is the vector of 
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direction cosines between zr and xcrYcrzcR* vector of direction 
cosines between z/^ and xcAycAZCA* Moreover, 0r,~0R, and 0 a are 3x3 matrices 
depending on oscillation of the axes x^yQZQ relative to XYZ, etc. Using equa- 
tions (10) and retaining only linear terms, the absolute velocities of typical 
points in the various substructures become 


Wc ' '^1 


9c '-2 9c 


«R “ 9 


R ^*2 9r 


«E = E, + Eg 

T ' T , T.T 


w« = An q. + A. q. 


( 11 ) 


'1 3A 

T . 


^2 ^A 
T . „T-,T 


* X^T^T’TX 

where = [wq^. I Qq ' DC^' ’ 9r = [woc! !c I DC < ®R^‘ ’ 9 e " ' §C ' DC i DE^ » 

and = [wQc! Qq j % ! §a 1 Da^^ configuration vectors for the substructures 


The system kinetic energy can be written in the form 




( 12 ) 


where 

^c ^ M ~c ~c ^ 1 9c ^11 §c 9c ^12 9c 9c ^22 9c 

Jmc 


T = — 
‘R 2 


1 


1 _T K- 


mr 


~R ~R '^"’r ^ 2 9r ^11 9r 9r *^12 9r 2 9r *^22 9r 


“ 2 


1 r -T 


1 iT 7 . 


m^ ~E '^'"e " 2 9e ^11 9e 9e ^12 9e 2 9e ^22 9e 


T _ 


03 


T = — 

‘a 2 


m„ 


'a ~a ^"^a 2 9a ^11 9a ^ 9a ^12 9a 2 9a ^22 9a 


in which 


c. 

- 

cl 

c,- 

dm^ , 

R.. = 

r{ 

Rn 

dm, 


i 

m J 
me 

1 

C 

13 J 

"■r 

1 

1 

E. 

_ 

eT 

Ei 

dm^ , 

A.. = ' 

''3 

m. 

A,- 

dm 

1C 

i j 

m^. ^ 

E 

U J 

1 

1 


(14 


Note that the square matrices C^j, IT-jj, E^j, and A'^-j have partitioned forms, 
with many of the off-diagonal suomatrices equal to zero. Introducing the n- 
dimensional configuration, vector for the entire system in the form g = 

tWoc ' §C ' ~R ' ~A IOC I De IDA^^’ where n is the number of degrees of freedom of 
the syitem^ the kinetic energy can be written in the general form 
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(15) 


T = ^q^Mq + q^Fq+lq^Kyq 

where M and Kj are symmetric matrices. Similarly, the potential energy for 
the entire system is 

V=lq’^K^q (16) 

where Ky is a symmetric matrix, and .the nonconservative virtual work has the 
form 

6W = Sq (17) 

where Q is the nonconservative generalized force vector. 

In general, the matrices M and F depend explicitly on time. However, un- 
der certain circumstances, such as when the substructures R and A are symmet- 
ric, the time dependence disappears. A helicopter with a symmetric rotor ro- 
tating relative to an airframe while in hover is an example, where the entire 
rotor is considered as a substructure. Another possibility is to consider 
each rotor blade as a separate substructure. In this case, a combination of 
substructures forms a symmetric rotor and M and F will once again be constant 
Tiatrices. 


Lagrange's equations can be written in the symbolic form 


^ Ik Ik n 

« sg ■ =*9 ' 8 


(18) 


i/here L = T - V is the system Lagrangian. Assuming that M and F are constant, 
introducing equations (15) and (16) into the Lagrangian L, and using equation 
’18), we obtain the Lagrange's equations of motion 

Mq + (f"^ - F)q + (K^ - K.^)q = Q (19) 


(here F^ - F is a skew symmetric matrix. Hence, equation (19) represents a 
ypical gyroscopic system. The natural frequencies and natural modes of the 
omplete structure and the closed-form solution of equation (19) can be ob- 
ained by the methods developed in references 12 and 13. The interest here 
s not so much in the response as in the dynamic characteristics of the sys- 
em, and in particular, the truncation effect on these characteristics. 


THE EIGENVALUE PROBLEM AND TRUNCATION IMPLICATIONS 


Introducing the 2n-dimensional state vector x(t) and the associated 2n- 
imensional force vector X(t) in the form 

x(t) = [ 9 ^(t) : 9 ^(t)f . X(t) = [9'^(t) 1 gV (20) 
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where 0 is the n-dimensional null vector, as well as the 2n x 2n matrices 


I = 


M 


0 . Kv 



(21) 


where 0 is the null matrix of order n, the n second-order differential equa- 
tions of motion, equation (19), can be replaced by the 2n first-order differ- 
ential equations in the state space x(t), where the equations have the eigen- 
value problem 


xlx + Gx = g (22) 

It is shown in reference 12 that the eigenvalue problem (22) can be reduced to 
the real symmetric form 

(o^Iy = Ky , i/lz = Kz (23) 

where K = gTI'Iq is a real symmetric matrix. The eigenvalue problem (23) is 
in terms of two real symmetric matrices and is known to possess real eigen- 
values. Assuming that I is positive definite, it follows that K is positive 
definite, so that the eigenvalues are not only real but also positive. More- 
over, the eigenvalues (r = l,2,...,n) have multiplicity two, so that to 
each belong the eigenvectors y^ and z^. Because I and K are positive defi- 
nite all the eigenvectors are independent. In fact, they are orthogonal with 
respect to the matrix I. 


Next, let us use the Cholesky decomposition and write I in the form I = 

LU , where L is a lower triangular matrix. Introducing the notation y^ = L'y^.. 
z^ = Mz^, (r = l,2,...,n), the eigenvalue problem (23) becomes 


y' 
r <r 




(24) 


where K' = L“^KL"T is a real symmetric positive definite matrix, in which 
l-t . (L-')T. 


Denoting by y an arbitrary 2n-vector, Rayleigh's quotient associated 
with the eigenvalue problem (24) can be written in the form (ref. 11) 

v^K'v 

R(v) = (25) 

y V 

Because K' is real and symmetric, it is well known that Rayleigh's quotient 
has a stationary value in the neighborhood of an eigenvalue. Note that the 
symmetric formulation (24) permits us to conclude that a stationarity prin- 
ciple exists also for gyroscopic systems. 
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Next, we wish to examine the truncation effect on the system characteris- 
tics. To this end, let us examine the eigenvalue problem Ay = xv, where A is 
a real symmetric matrix of order N, and assume that the eiginvalues of A are 
ordered so that x^ £ ^2 — ••• — ^N- let us form the matrix B by deleting 

the last row and column of A and write the eigenvalue problem By = yy> where 

the eigenvalues Yj (j = 1>2 N-1) are ordered so that y] 1. Yo £ ••• £ Ym_T • 

The question arises as to how the eigenvalues yj relate to the eigenvalues x^-. 
To this end, one can use the Courant's maximum-rninimum theorem (ref. 14) and 
prove that 


^1 i^l 


^2 - 


< X 


N-1 


^N-1 - ^N 


(26) 


We shall refer to inequalities (26) as the inclusion principle . 


Now, let us return to the truncation problem. The 2n x 2n matrix K' was 
obtained as the result of representing the spinning structure by an n-degree- 
of-freedom system. Note that the rotational coordinates are also included in 
these degrees of freedom. This representation is tantamount to the imposition 
of a given number of constraints on the original structure. For example, the 
first of equations (10) can be written in the form 


He “ X ?Ci 


i=l 


Ci ’^Ci 


(27) 


so that the constraints imposed on the system are nc,nc+l “ ^C,nc+2 ~ ~ 

Truncating the series (27) by assuming that nc,nc " obtain a matrix K" 

obtained from K' by deleting two rows and the corresponding two columns. If 
the eigenvalues of K' are such that £ ^ 2 ^ £ .a* £ eigen- 

values Bp of K" are such that B^^ £ B 2 ^ £ ••• £ ^n-1 ’ 

9 9 9 9 9 9 9 

0 )^ £ B^ £ t 02 £ B 2 £ . • . £ (^n -1 - ^n-1 - “n 


Note that the fact that the eigenvalues of K' and K" have multiplicity two is 
automatically taken into account in inequalities (28). On the other hand, by 
relaxing one constraint, i.e., by adding one term to the series (27), we ob- 
tain a (2n + 2) X (2n + 2) matrix K'" which is obtained by adding two rows 
and columns to K'. The eigenvalues of K'" are such that 


2 2 2 2 2 2 2 
“■j £ ^^] £ “2 — ■ — °^n — ‘*’n — °n+l 


(29) 


The above developments permit us to conclude that the system estimated natural 
frequencies tend to decrease monotonically with each additional degree of 
freedom. At the same time there is a new frequency added which is higher 
than any of the previous ones. 

The question remains as to how to select the admissible functions or ad- 
missible vectors. The first thing that comes to mind is to take them as the 
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eigenfunctions and eigenvectors of the various substructures. In many cases, 
the solution of the eigenvalue problem for a substructure can be quite a 
task in itself, so that in such cases one may wish to use deformation patterns 
only approximating the actual modes. This can be regarded as imposing addi- 
tional constraints on the system, which tends to raise the natural frequencies 
of the system, but this may be considered as a viable alternative, particular- 
ly when the validity of the solution of the eigenvalue problem is question- 
able. Experience with the Rayleigh-Ritz approach shows that the system natu- 
ral frequencies are not very sensitive to the admissible functions used, which 
can be traced to the stationarity principle. But a stationarity principle 
exists also for discrete systems, so that the same conclusion can be extended 
to admissible vectors. 

The truncation by substructures has a clear advantage over truncation 
of the structure as a whole. The reason is that it permits a more rational 
judgement based on the substructure properties, such as the mass and stiffness 
distributions. Generally one is interested in only a limited number of lower 
modes of the complete structure. Hence, a very stiff and light substructure 
is likely to have less effect on the modes of the complete structure than a 
flexible heavy substructure. Hence, one can truncate the first more severely 
than the second. Some ideas for truncation can be obtained by estimating the 
natural frequencies of the substructures. This by no means implies that one 
need solve the eigenvalue problem for the substructures exactly. Indeed, 
using a Rayleigh-Ritz procedure for continuous or discrete systems, in con- 
junction with a preselected set of admissible functions or admissible vectors, 
it is possible to obtain a reasonable estimate of the lower frequencies of 
each substructure. Note that the Rayleigh-Ritz method can be used to produce 
and solve an eigenvalue problem of considerably lower dimension than that of 
the full eigenvalue problem for the substructure. Tho estimated lower natural 
frequencies of the substructure, when compared to those of other substructures, 
can be used merely as a guide for truncation purposes. In fact, the eigen- 
vectors serve no useful purpose and need not be calculated, as the same ad- 
missible functions or vectors can be used to represent the substructure in 
the generation of the eigenvalue problem for the complete assembled struc- 
ture. This conclusion is based on results shown in reference 11. 

If the dimension of the eigenvalue problem for the complete assembled 
structure is still too large, and the higher modes are not really necessary, 
then one can solve only for a given number of lower modes by using such tech- 
niques as subspace iteration. 


CONCLUDING REMARKS 


A procedure has been shown whereby the equations of motion for large 
structures with rotating flexible components can be derived by the Lagrangian 
approach. A fundamental consideration in the derivation of Lagrange's equa- 
tions is the superposition of substructure motions by means of an orderly 
kinematical procedure, which automatically eliminates the problem of con- 
straints. Using a Rayleigh-Ritz approach, it is shown that each continuous 
or discrete flexible substructure can be simulated by a finite number of ad- 
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missible functions or admissible vectors and exact substructure modes are 
not really necessary. This conclusion is based on a stationarity principle 
for rotating structures developed recently by the first author (ref. 11). Fi- 
nally, the substructure synthesis approach provides a rational basis for trun- 
cating the number of degrees of freedom both of each individual substructure 
and of the assembled substructure. 
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